function [ x ] = elim_gauss( A, b )
%resuelve Ax=b usando eliminacion de Gauss sin pivoteo
%A es matriz cuadrada y b vector vertical
    extendida = [A b];
%triangulacion
    tam_A = size(A,1);
    for j=1:(tam_A-1)
        for i=(j+1):tam_A
            extendida(i,:) = extendida(i,:) - extendida(i,j) * extendida(j,:) / extendida(j,j);
        end
    end
%sustitucion hacia atras
    x = zeros(tam_A,1);
    for i=tam_A:-1:1
        sumatoria = 0;
        for j=tam_A:-1:(i+1)
            sumatoria = sumatoria + extendida(i,j)*x(j);
        end
        x(i) = ( extendida(i,tam_A+1) - sumatoria ) / extendida(i,i);
    end
end